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1. INTRODUCTION 

The influence of the earth's magnetic field on the distribution of 
ionization in the ionosphere is well-known. For example, the equatorial 
anomaly is roughly symmetrical with respect to the dipole equator rather 
than the geographic. Furthermore, radio occultation measurements of the 
Pioneer' and Voyager space probes indicate the possibility of the 
existence of equatorial anomaly in the ionospheres of Jupiter (Mahajan, 
1981) and Saturn (Mahajan, et al . , 1985). 

The formation of equatorial anomaly in the Jovian ionosphere under 
model electric fields was investigated by Tan (1986). Tan (1990) further 
studied the differences between the topside ionospheric profiles in the 
Eastern and Western Sectors of Jupiter due to the offset effect of its 
eccentric dipole. 

In this study we investigate the offset effect of Saturn's dipole 
on its ionosphere. The magnetic field of Saturn is primarily that of a 
dipole closely al igned to the rotational axis but displaced northward 
from the center by a distance approximately equal to 0.05 R $ , R s being 
the reference radius of Saturn (cf. Connerney, et al . , 1984). This offset 
effect would manifest itself most prominently between the ionospheric 
profiles in the Northern and Southers hemispheres of Saturn. 



2. SATURN'S DIPOLE MODEL 

The magnetic field of Saturn can be approximated by that of a 
dipole coaxial with the rotational axis and displaced from the center 
northward by a distance of a = 0.05 R s (Connerney, et al . , 1984). The 
following systems of coordinates are defined in a meridional plane: 

(1) Planetocentric cartesian (x 1 , z‘), 

(2) Planetocentric polar (r 1 , 8'), 

(3) Dipole cartesian (x, z), and 

(4) Dipole polar (r, 9) coordinates. 

The transformation from (r, 0) to (r‘ , 9') is effected in the following 
stages: (r, 0) -A-* (x, y) -B+ (x', z') -C- (r‘ , 9'), where the 
transformations A, B and C are given by (vide Fig . 1 ) 

A: x = r sin9 , 

z = r cos9 , 

B: x' = x , 

z' = z + a , 

C: r' = (x' 2 + z' 2 )* . 

0' = tan“^(x'/z') , 

Substituting Eq. (2.1) - (2.4) in (2.5) and (2.6), we get 

r' = (r 2 + 2ar cosQ + a 2 )* , 

^ -1 , sin0 , 

and 9 = tan ( cos6 + ■ 

The inverse transformations from (r 1 , 9') to (r, 0) are obtained by 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 
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interchanging the primed and unprimed quantities and reversing the sign 
of a. 


r = (r 1 2 - 2ar' cos0' + a 2 )* 


(2.9) 


and 


tan 


/ sine 1 
''cose 1 - a/r 


r) 


( 2 . 10 ) 


For a dipole field line of definite L-vaule, the coordinates of 

the foot of the field line in the northern hemisphere (r Q , e Q ; r' Q , Q' q ) 

can be cal culated given the radial distance r' Q = R s + h o = R ( say )’ where 

h is the altitude of the base of the ionosphere. For a diffusive 
o 

equilibrium model, h Q is approximately the height of the maximum electron 
density. For Saturn, h Q = 2,500 km (from Mahajan, 1987) and R $ = 58,300 
km, whence R = 60,800 km. 

The foot of the field line is the intersection of the circle r‘ = R 
and the dipole field line, whose equation is 


r = L R^ sin 2 e . (2.11 ) 

Applying Eq. (2.11) in (2.7), we get 

L 2 R<. 2 sin 4 e Q + 2aL R s sin 2 e Q cose o + a 2 - R 2 = 0. (2.12) 

One can solve for y = sine o numerically by the Newton-Raphson scheme 
(cf. Stark, 1970) 


f(Yi) 

y i+l = y i " "FTTjT ’ 

where f(y^) = L 2 R^ 2 y^ 4 + 2aL R^ y^ 2 0 - Yj )* + a R » 

, ? , yi (2 - 3y, 2 ) 

and flyp - 4 L 2 R $ 2 y, 3 + 2aL R $ - - • 


(2.13) 

(2.14) 


(2.15) 
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3. THEORY 

The topside ionospheric plasma can be approximated by a static 
diffusive equilibrium model V g = 0, where V $ is the velocity of the ions 
resolved along the field line (cf. Bauer, 1969). This is even more true 
for the ionospheres of the major planets, where the production and loss 
rates are extremely slow. V g is prescribed by the momentum equation 

V s = / L m 9 s + mf s - -jr ( NkT) ] , (3.1) 


where N = number density of electrons/ions, 
m = mass of ions (H + ), 
v = ion-neutral collision frequency, 

T = Temperature, 
k = Boltzmann's constant, 
s = arc length along dipole field line, 
g $ = gravitational acceleration resolved along s, 
and f s = centrifugal acceleration resolved along s. 

Under the condition of static diffusive equilibrium, Eq. (3.1) integrates 
easily from the foot of the field line in the northern hemisphere to the 
field point 


N 



.s m 

a f oW? 

e 


(g s + f s ) ds 


(3.2) 


The differential arc length is given by 


ds = L R $ (3 cos 2 0 + 1)* sine de . 

Taking into account the radial dependence of g, we have 


(3.3) 
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9 0 R 

..2 


sinl 1 


(3.4) 


where I* is the true inclination of the field line and g Q the acceleration 
of gravity at the foot of the field line in the northern hemisphere. 

From Fig. 2, we have 

I' = I + 9' - 9 


= tan -1 (2 cot 0 ) + tan 1 ( 


sin9 ) - tan" 1 (tan 0 ) . (3.5) 


cos 0 + a/r 


Also from Fig. 2, 


f = n 2 r 1 sin 0 ' cos(I' - 0 ' ) . 

s 


(3.6) 


Since r* sin 0 ' * r sin 0 , I' - 0 ' = I - 9 and 


cos( I - 0 ) = 


3 cos9 sin 6 
(3 cos 2 0 + l ) 4 


(3.7) 


we get, by substituting from ( 2 . 11 ) 

4 

3 sin 9 cos 0 

f s = Q L R S (3 cos 2 0 + 1)* * 


(3.8) 
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4. INPUT DATA 

In the diffusive equilibrium model, it is appropriate to take the 
electron density at the base as the observed peak electron density. In 
the present study, this is assumed to be 10 electrons/cm (from Mahajan, 
1987). The altitude of the base is taken as the observed height of the 
peak, which is roughly 2,500 km (Mahajan, 1987). 

The electrons, ions and neutrals are assumed to be in thermal 
equilibrium (Henry and McElroy, 1969; Nagy, et al., 1976). The temperature 
profile is taken from the measurements of Festou and Atreya (1982), who 
found an isothermal region above 1450 km of constant temperature of 800 K. 

The reference radius of Saturn as well as the angular velocity of 
rotation are taken from Mullin (1984). 
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5. RESULTS AND DISCUSSIONS 

The diffusive equilibrium model is most appropriate for equinox 
conditions when both seasonal and ring effects are absent. This justifies 
the assumption that the base electron densities in the two hemispheres 
are equal . 

Figures 3 through 5 illustrate the results of our calculations. 

In Fig. 3, contour plots of the electron density in a meridional plane 
in the two hemispheres are shown together with the dipole field lines. 

In general, the electron densities are higher in the northern hemisphere 
than in the southern. This is a direct consequence of the dipole field 
line alignments in the two hemispheres. In the northern hemisphere, 
ionization is transported along the field lines to relatively lower 
altitudes than in the southern hemisphere, thus producing enhanced 
electron densities in the topside ionosphere. 

Also shown in Fig. 3 are the loci of points at which the components 
of gravitational and centrifugal forces along the field lines are equal 
and opposite. In the region of centrifugal dominance, the electron 
density slowly increases with distance, thus creating a region of 
enhanced electron density away from the planet. As a result, troughs of 
electron density are found at midlatitude heights (X = 25° to 35° in the 
northern hemisphere and X = 25° to 30° in the southern hemisphere) rather 
than above the equator. Above the equator, the electron density profile 
possesses a minimum which is the saddle point in the electron density 
contour. 

Figures 4 and 5 show the vertical profiles of electron density in 
the two hemispheres at a midlatitude (X = 30°) abd a low-latitude (X - 15 ) 
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location. As stated earlier, the electron densities were generally 
greater in the northern hemisphere from the same base electron density. 
Consequently, the scale heights are also greater in the northern 
hemisphere. 

The electron densities at the midlatitude locations were smaller 
than those at the low-latitude locations. In fact, the troughs in the 
electron density contour plots are located around roughly 30° latitudes 
in either hemisphere. At the low latitude locations, the electron density 
profiles undergo minima and actually Increase with altitude thereafter. 
This would be true so long as the assumption of perfect corotation of 
plasma remained valid. 

The calculated scale heights in the topside ionosphere for the 
latitudes of 15° and 30° in both hemispheres are given in Table I. The 
difference in the scale heights between the two hemispheres of around 
200 km is significant. This difference should reveal itself if sufficient 
observed data become available. Also shown in the Table are the observed 
scale heights deduced from Voyager 2 measurements (from Mahajan, 1987) 
which correspond to near equinox conditions. The calculated scale heights 
were significantly greater than those observed, which are about 1000 km. 
The scale height calculated from model computations of Atreya and Waite 
(1981) is about 1050 km. 
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TABLE 

I. CALCULATED AND OBSERVED SCALE 

HEIGHTS 

CALCULATED SCALE HEIGHTS 

LATITUDE 

HEMISPHERE 

SCALE HEIGHT 

15° 

N 

2,000 km 

15° 

S 

1 , 800 km 

30° 

N 

1,700 km 

O 

O 

CO 

S 

1,500 km 

OBSERVED SCALE HEIGHTS* 

LATITUDE 

HEMISPHERE 

SCALE HEIGHT 

36° 

N 

1,000 km 

31° 

S 

1,000 km 


* From Mahajan (1987) 


X 



Fig.l. Planetocentric Cartesian (x',z'), Planetocentric Polar (r',0 1 ), Dipole Cartesian (x,z) and 
Dipole Polar (r,0) Coordinates. 


£L 2 r'sin 6 1 cos (1-09 



Fig. 2. Components of gravitational and centrifugal accelerations resolved along the dipole field line. 



Fig. 3. Contour Plot of electron density in the meridional planes of Northern and Southern Hemispheres 
for equinox conditions. The solid lines are dipole field lines whereas the broken lines are 
contours of constant electron density. The dotted lines are loci of points at which the 
components of gravitational and centrifugal forces along the field lines are equal and opposite. 
The cross mark represents the saddle point in the electron density profile. 
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Saturn's Plasmasphere in Offset Dipole Model* 


1.1991 Spring Meeting 


A Tan (Department of Physics, Alabama A&M University, 
Normal, AL 35762; 205-851-5313) 

The magnetic field of Saturn can be approximated by 
that of a dipole coaxial with the rotational axis and 
displaced from the center northward by a distance of 
0.05 Rg, Rg being the reference radius of Saturn. A 
diffusive equilibrium model of Saturn's ionosphere for 
equinox conditions is constructed to study the offset 
effect of this dipole. Because of the field line 
alignments, the electron densities are generally 
greater in the northern hemisphere than in the 
southern. The plasma scale heights are also greater in 
the northern hemisphere. There is plasma pileup at 
higher heights above the equator, owing to the 
centrifugal dispersion of ionization there. 
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